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Abstract We present a study of the generalized parton 
distributions (GPDs) for the quarks in a proton in both 
momentum and position spaces using the light-front 
wave functions (LFWFs) of a quark-diquark model for 
the nucleon predicted by the soft-wall model of AdS/ 
QCD. The results are compared with the soft-wall AdS/ 
QCD model of proton GPDs for zero skewness. We also 
calculate the GPDs for nonzero skewness. We observe 
that the GPDs have a diffraction pattern in longitudi¬ 
nal position space, as seen before in other models. Then 
we present a comparative study of the nucleon charge 
and anomalous magnetization densities in the trans¬ 
verse plane. Flavor decompositions of the form factors 
and transverse densities are also discussed. 


1 Introduction 

Hadronic structure and their properties being nonper- 
turbative in nature are always very difficult to evaluate 
from QCD first principle and there have been numer¬ 
ous attempts to gain insight into hadrons by study¬ 
ing QCD inspired models. The quark-diquark model, 
where a nucleon is considered to be a bound state of 
a single quark and a scalar or vector diquark state, is 
proven to reproduce many interesting properties of nu¬ 
cleons and has been extensively used to investigate the 
proton structure. Recently, a light front quark-diquark 
model for the nucleons has been proposed in Ref. pQ, 
where the light front wave functions are modeled by 
the wave functions obtained from a soft-wall model in 
light front AdS/QCD. The light front wave functions 
(LFWF) are derived by matching the electromagnetic 
form factors of hadrons in the light front QCD and soft- 


wall model of AdS/QCD. The model is consistent with 
Drell-Yan-West relation relating the high Q 2 behavior 
of the nucleon form factors and the large x behavior 
of the structure functions. Recently, LFWF of baryon 
and the light baryon spectrum have been described by 
extending the superconformal quantum mechanics to 
the light front and embedding it in AdS space [2]. The 
LFWF for the rho meson in AdS/QCD has been suc¬ 
cessfully applied to predict the diffractive rho meson 
electroproductionpj]. 

In this paper, we study the proton structure and 
evaluate the Generalized Parton Distributions(GPDs), 
transverse charge and magnetization densities in the 
light front quark-diquark model. Contrary to ordinary 
parton distribution function, GPDs are functions of three 
variables, namely, longitudinal momentum faction x of 
the quark or gluon, square of the total momentum trans¬ 
ferred (f) and the skewness £, which represents the lon¬ 
gitudinal momentum transferred in the process and pro¬ 
vide interesting information about the spin and orbital 
angular momentum of the constituents, as well as the 
spatial structure, of the nucleons (see [¥j for reviews on 
GPDs). The GPDs appear in the exclusive processes 
like deeply virtual Compton scattering (DVCS) or vec¬ 
tor meson productions and they reduce to the ordinary 
parton distributions in the forward limit. Their first 
moments are related to the form factors and the sec¬ 
ond moment of the of the sum of the GPDs are related 
to the angular momentum by a sum rule proposed by 
Jig]. Being off-forward matrix elements, the GPDs have 
no probabilistic interpretation. But for zero skewness, 
the Fourier transforms of the GPDs with respect to the 
transverse momentum transfer (Z\j_) give the impact 
parameter dependent GPDs which satisfy the positivity 
condition and can be interpreted as distribution func¬ 
tions [Q. The transverse impact parameter dependent 
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GPDs provide us with the information about partonic 
distributions in the impact parameter or the transverse 
position space for a given longitudinal momentum (a;). 

The impact parameter b± gives the separation of the 
struck quark from the center of momentum. In paral¬ 
lel to the efforts to understand the GPDs by theoreti¬ 
cal modeling, different experiments are also measuring 
deeply virtual Compton scattering and deeply virtual 
meson production to gain insight and experimentally 
constrain the GPDs[7j. 

We evaluate the proton GPDs for both zero and 
nonzero skewness and compare with the results in a 
soft-wall AdS/QCD model (8|(for hard-wall and soft- 
wall AdS/QCD models of hadrons, see |9llT0]'l . For 
zero skewness, the GPDs are investigated in the im¬ 
pact parameter or transverse position space. The LF 
diquark results for GPD H(x,b±_) for u-quark is almost 
the same as AdS/QCD results whereas there is a lit¬ 
tle difference for d-quark. But the LF diquark model 
results for E(x,b±) for both u and d quarks are differ¬ 
ent from AdS/QCD results. For nonzero skewness, the 
GPDs in longitudinal impact parameter space show a 
diffraction pattern. It is interesting to note that sim¬ 
ilar diffraction patterns were observed in simple QED 
model for DVCS amplitude jTl] and GPDs[T2j and in a 
phenomenological model of proton GPDsjTS]. 

Electric charge and magnetization densities in the 
transverse plane also provide insights into the structure 
of nucleons. The charge and magnetization densities in 
the transverse plane are defined as the Fourier trans¬ 
form of the electromagnetic form factors. The form fac¬ 
tors involve initial and final states with different mo¬ 
menta and the three dimensional Fourier transforms 
cannot be interpreted as densities whereas the trans¬ 
verse densities(i.e., Fourier transformed only for trans¬ 
verse momenta) defined at fixed light front time are free 
from this difficulty and have proper density interpretation [T¥l 
US]. We calculate the transverse charge and anoma¬ 
lous magnetization densities for both proton and neu¬ 
tron in the light-front diquark model and compare with 
the two different global parameterizations proposed by 
Kelly [H] and Bradford et al [17]. We present results 
for both unpolarized and transversely polarized nucle¬ 
ons. We also present a comparison with the AdS/QCD 
results citeCM3 for the transverse charge and magneti¬ 
zation densities. 

The paper is organized as follows. In Section [2] we 
give a brief introductions about the nucleon LFWFs of 
quark-diquark model as well as the electromagnetic fla¬ 
vor form factors. We show the results for proton GPDs 
of u and d quarks in momentum space in Section [3] 

Then we discuss the GPDs in the transverse as well 
as the longitudinal impact parameter space in Sections 


urn and We present the results of the charge and 
anomalous magnetization densities in the transverse plane 
in section [4] Finally we provide a brief summary and 
conclusions in Section [5] For GPDs with nonzero skew¬ 
ness, we present a comparison of the quark-diquark 
model results with a Double Distribution(DD) model 
in the appendix. 


2 Light-front quark-diquark model for the 
nucleon 

In quark-scalar diquark model, the three valence quarks 
of nucleon are considered as an effectively composite 
system composed of a fermion and a neutral scalar 
bound state of diquark based on one loop quantum fluc¬ 
tuations. In the light-cone formalism for a spin A com¬ 
posite system the Dirac and Pauli form factors F\{q 2 ) 
and F^q 2 ) are identified to the hclicity-conserving and 
helicity-flip matrix elements of the J + current [19] 


(P 
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+ 0,t 
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( 2 ) 


here M n is the nucleon mass. Writing proton as a two 
particle bound state of a quark and a scalar diquark 
in the light front quark-diquark model, the Dirac and 
Pauli form factors for the quarks can be written in the 
light-front representation as 
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,(4) 


s.j_ — kj_ + (l—x)qj_. ip\” q (x, kj_) are the LFWFs 
with specific nucleon helicities A n = ± and for the 
struck quark X q = ±, where plus and minus corre¬ 
spond to and — / respectively. We consider the 
frame where q = (0, 0, qj_), thus Q 2 = —q 2 = q/_. 

We adopt the generic ansatz for the quark-diquark 
model of the valence Fock state of the nucleon LFWFs 
at an initial scale fj ,o = 313 MeV as proposed in PQ : 


4>+ q {x, kj.) = ^Obkx), 
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k 1 + ik 2 
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Ip + Jx, k ± ) = 


k 1 — ik 2 


<?l 2) 0,k ±), 


xM n *® 

^Z,(a;,kj.) = ^ 1} (o:,kj.), (5) 

where ^^(x, kj_) and (fiq 2 \x, kj_) are the wave func¬ 
tions predicted by soft-wall AdS/QCD [51], modified by 
introducing the tunable parameters and bq' 1 for 
quark q [I]: 

Av V J- 3/ 


exp 


log(l/x) 


( 6 ) 


2re 2 (1 — x) 2 

The parameters are tuned to fit the electromagnetic 
properties of the nucleons. Following the convention of 
PQ, we fix the normalizations of the Dirac and Pauli 
form factors as 


F?(Q 2 ) = n, 


nm „ nm 


F% iQ 2 ) = «, 


2^ y-^-jKo) > 


(7) 


9 I 9 (0) ’ 

so that F/(0) = n 9 and F|(0) = where n„ = 2 , rid = 
1 and the anomalous magnetic moments for the u and d 
quarks are k u = 1.673 and re^ = —2.033. The advantage 
of the modified formulae in Eq. 0 is that, irrespective 
of the values of the parameters, the normalization con¬ 
ditions for the form factors are automatically satisfied. 
The structure integrals, J 9 (Q 2 ) have the form as 

IiiQ 2 ) = [ dxx 2 ^ (1 - x) 1+2b ^ R q (x, Q 2 ) 
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with 


R q (x, Q 2 ) = 1 + cr 2 (x) 
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( 10 ) 


It is straightforward to write down the flavor decompo¬ 
sitions of the Dirac and Pauli form factors of nucleon 


as 

F P^) =euF uW +edF d(u)^ 


(* = 1 , 2 ) 


( 11 ) 


where e u and are the charges of u and d quarks in 
units of positron charge(e). On top of the AdS/QCD 
scale parameter re, the wave functions involve four more 
parameters and bq 1 (with i = 1,2) for each quark. 


Parameters 

U 

d 

oW 

0.035 

0.20 

fed) 

0.080 

1.00 

ad) 

0.75 

1.25 

fed) 

-0.60 

-0.20 

iv d) 

29.180 

33.918 

ivd) 

1.459 

1.413 

Table 1 The parameters in the model for re = 0.4066 GeV 

Quantity 

Our results 

Measured data [23] 

r P 
' E 

0.7861 fm 

0.877 ± 0.005 fm 

r P 
' M 

0.7719 fm 

0.777 ± 0.016 fm 

(r%) n 

-0.085 fm 2 

-0.1161 ±0.0022 fm 2 

r n 

' M 

0.7596 fm 

0.862±°;°°® fm 

Table 2 electromagnetic radii of the nucleons 

In Ref.pQ, re is 

taken to be 0.35 GeV and the parame- 

ters are evaluated to fit the electromagnetic properties 
of the nucleon. But the results for the form factors pre- 

sented in that 

paper are 

not converged with respect 


to the lower limit of x integrations in Eqs.flS] and 0. 
The comparisons with experimental data presented in 
several plots in Ref.fT] are true only for an unrealisti¬ 
cally large value of lower limit for the x integrations 
which drastically change when integrated from x —> 0. 
So, when proper limits in the x-integrations are taken, 
the parameters presented in [I] cannot reproduce the 
data. In this work, we use a different scale parameter 
re = 0.4066 GeV which was obtained by fitting the nu¬ 
cleon form factors in AdS/QCD soft-wall model [8][55]. 
Here, we show that we can reproduce the nucleon form 
factors with the new parameters and IVW listed 

in Table [T] The results are stable and converged under 
the integration over x. The new parameters reproduce 
the experimental data quite accurately for a wide range 
of Q 2 values. 

The Dirac and Pauli form factors of both u and d 
quarks are shown in Fig [I] The form factors F\ 9 and F% 
for both u and d quarks in light-front quark-diquark 
model for the scale parameter re = 0.4066 GeV and the 
parameters defined in Table [I] are in excellent agree¬ 
ment with the data. For Ff, we can see a clear improve¬ 
ment in the quark-diquark model over the AdS/QCD 
model. It is important to note that other models fail 
to reproduce the form factors data for d quark [Mj. In 
Fig[21 we have shown the fit of light-front quark-diquark 
model results with experimental data of proton form 
factors. We get excellent agreement with the data. In 
the same plots, we also show comparisons of the light- 
front quark-diquark model and the soft-wall AdS/QCD 
model with the same value of re [52|. The results of 
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Fig. 1 (Color online) Plots of flavor dependent form factors for u and d quarks. The experimental data are taken from 1271 
1281 . The red dashed lines represent the soft-wall AdS/QCD model 1221 . 




Fig. 2 (Color online) Light-front quark-diquark model results are fitted with the experimental data. The plots show the ratio 
of Pauli and Dirac form factors for the proton, (a) the ratio is multiplied by Q 2 = —q 2 = — t, (b) the ratio is divided by k p . The 
experimental data are taken from Refs. I29II30|131II32II33| . The red dashed lines represent the soft-wall AdS/QCD model [22). 


the light-front quark-diquark model agree with the data 
better than AdS/QCD, specially at large Q 2 values we 
achieve substantial improvement. The Sach form factor 
G e (Q 2 ) for the neutron is shown in FigGD Again, our 
results agree with the experimental data much better 
than the AdS/QCD results. The fitted results for the 
electromagnetic radii of the nucleons are listed in Table 


[21 The standard formulae for the electromagnetic radii 
of nucleon used here are given below: 
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(r 2 M ) N = 
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Fig. 3 (Color Online) The Sach form factro G^(Q 2 ) for the neutron. The experimental data are taken from Refs. 134113511361 
. The red dashed line represents the soft-wall AdS/QCD model prediction. 


where N stands for nucleon(iV = p/n) and the Sachs 
form factors are defined as 

g n e(Q 2 ) = F xiQ 2 ) - 7Sr*f(Q 2 ), 


4 M 2 N 

N ( 


GZ(Q 2 )=F?(Q') + F?(Q'). 


(14) 

(15) 


3 Generalized parton distributions 

Using the overlap formalism of light front wave func¬ 
tions, we evaluate the GPDs in light front quark-diquark 
model. We consider the DGLAP domain , i.e., £ < x < 
1 where £ is the skewness and x is the light front longi¬ 
tudinal momentum fraction carried by the struck quark. 
This domain corresponds to the situation where one re¬ 
moves a quark from the initial proton with light-front 
momentum fraction x = -j=q- and the transverse mo¬ 
mentum ku and re-insert it into the final state of the 
proton with longitudinal momentum fraction x — £ and 
transverse momentum — qj_- The contributions to 
the GPDs for 0 < x < C, come from the particle num¬ 
ber changing interactions and cannot be studied in this 
model. The kinematical domain for GPDs studied here 
is thus restricted to ( < x < 1 where only diagonal 
overlaps(2-particle state —> 2-particle state) contribute. 
The GPDs H and E are defined through the matrix el¬ 
ement of the bilocal vector current on the light-front: 

f dy ~ A* 

./ 87T 


■ /2 <P',A'| riOh+^y) \P, A) 


1 


2P+ 






E(x,C,t) 


2M, 


-<J +a q a 


U(P, A). 


(16) 


The proton state | P, A) is written in two particle Fock 
states with one fermion and a scalar boson in the light 
front quark-diquark model. Using the relations 


1 t/(P', X'h + U(P, A) = § x x , 


1-1 


2P+ 


4(i 


S\,X' + 


1 Aq 1 + iq 2 
2 M n 


(17) 


where P = (P + P')/2 and A(A') = ±| is the ini- 
tial(final) proton spin, we have the following expressions 
for the GPDs in terms of the LFWFs in the quark- 
diquark model 
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(c) 



(b) 


(d) 




Fig. 4 (Color online) Plots of (a) H u (x,t) vs x (b) the same as in (a) but for d quark (c) E u (x,t) vs x (d) the same as in (c) 
but for d quark for fixed values of —t. The black lines represent the light-front quark-diquark model and the red lines are for 
AdS/QCD model 0. 


Substituting the LFWFs (Eq.flSJ) in Eos. lfTHl) and (fTTfl) 
and integrating over k^, we get the following expres¬ 
sions for GPDs 

,t) 
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where Q 2 = —1( 1 — £) — Tf 2 £ 2 and -<4 is a function of x 
and x': 

log(l/x) log(l/x') 


A = A(x, x’) = 


(25) 


2k 2 (1 — x) 2 2k 2 (1 — x') 2 

/■f(O) and /f (0) are the integrals defined in Eqs. (JHJ) and 
@ for Q 2 = 0. The GPDs are normalized as 

f dxH q (x : 0,0) = n q , 
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where n q denotes the number of u or d valence quarks 
in the proton and the quark anomalous magnetic mo¬ 
ment is denoted by n q . According to the polynomiality 
condition, the n-th Mellin moment of a GPD should 
be a polynomial with highest power at t —> 0(25; 
[20j. Since the moments require the GPDs for all val¬ 
ues of —1 < x < 1, it is not possible to confirm the 
polynomiality condition in this model as the GPDs are 
evaluated only for 0 < £ < x. But we have numeri¬ 
cally checked for n = 1,2,3, and 4 that the moments 

x n ~ 1 H(x, £, t) show the behavior consistent with the 
polynomiality condition in the limit t —> 0. 

The GPDs for zero skewness(£ = 0) in light-front 
quark-diquark model are compared with the AdS/QCD 
results [1] in FigQ] In Figs. Oja) andQJb), we show the 
GPD H ( x , t) as a function of x for different values of —t 
for u and d quarks. Similar plots of E{x, t) for u and d 
quarks are shown in Figs0Jc) and0Jd). The overall na¬ 
ture of both the models is same for u quark while there 
are some disagreements in the GPD H(x,t) for the d 
quark. Since the d-quark form factors are not well de¬ 
scribed in AdS/QCD, this disagreements are expected. 
The GPD H(x,t ) falls off faster as x increases for d 
quark compare to u quark in both the model. Unlike 
H(x,t), the fall-off of the GPD E{x,t) at large x is 
similar for both u and d quarks with increasing x. 

In Fig[5l we show the skewness dependent GPDs as 
a function of x and t for a fixed ( = 0.15. The overall 
behaviors of the GPDs with nonzero £ are similar to 
the zero skewness GPDs. The same GPDs for a fixed 
x = 0.8 are shown as a function of ( for different values 
of —t in Fig. [G] In FigJTfa) andjTJb) we have plotted the 
GPD H (a;, £, t ) as a function of x for u and d quarks for 
different values of ( with fixed value of —t = 0.8 GeV 2 . 
The similar plots of E{x, /, t) for u and d quarks are 
shown in FigJTJc) and Eld). In FigEl the peaks of all 
the distributions move to higher x as £ increases and the 
amplitudes of the distributions increase with increasing 
£ for a fixed value of —t. Due to the factor of y/1 — £ in 
the denominator of the the GPD id(x, £,f) in Ea. E^l) . 
the increase in the magnitude with increasing £ is more 
in H(x,£,t) than in E{x,C,,t). 


3.1 GPDs in transverse impact parameter space 

GPDs in transverse impact parameter space are defined 
as [6ll43|: 

H{x, £.b ± ) = I d 2 ^e-^H { x M 
E(x, £, bx) = J d 2 q ± e-^ h ^E(x, £, t). (27) 


Here, bx is the transverse impact parameter. For zero 
skewness, bx gives a measure of the transverse distance 
between the struck parton and the center of momentum 
of the hadron. bj_ satisfies the condition J2i x ib±i = 0, 
where the sum is over the number of partons. The rela¬ 
tive distance between the struck parton and the center 
of momentum of the spectator system is given by , 
which provides us an estimate of the size of the bound 
state 03. However, the exact estimation of the nuclear 
size is not possible as the spatial extension of the spec¬ 
tator system is not available from the GPDs. In the 
DGLAP domain x > Q the impact parameter bx im¬ 
plies the location where the quark is pulled out and 
pushed back to the nucleon. In the ERBL region x < (, 
bx gives the transverse location of the quark-antiquark 
pair inside the nucleon. 

In Fig 1H1 we compare the transverse impact param¬ 
eter dependent proton GPD H(x, bx) for zero skewness 
in the light-front quark-diquark model and in the soft- 
wall AdS/QCD. Unlike H d (x,h±), the diquark model 
results for H u (x, bx) are in good agreement with AdS/ 
QCD. The GPD H d (x,h±) fall off slowly for large x 
in AdS/QCD compared to the diquark model while the 
fall-off of H u (x, bx) in both models is same. The reason 
of the disagreement in H d {x, bx) is that the AdS/QCD 
model is unable to reproduce F d to match with exper¬ 
imental data whereas the form factor in the diquark 
model agrees well with the data (see FigQJb)). The 
overall shapes of the curves in FigEDa) and (c) are due 
to the fact that the two particle LFWFs are effectively 
functions of “x-weighted transverse variable” [21] z = 
\/x(l — x)|&x| which is true for both the AdS/QCD 
and the diquark models. Both the models lack the sym¬ 
metry about x —> (1—x), the asymmetry is more promi¬ 
nent for d quark, in the the diquark model due to the 
parameters aq d ^ bq' 1 in the diquark wave functions. 
In Fig. [9] we have compared the two models for the 
proton GPD E(x,h±). The GPD E(x,h±) in impact 
parameter space in the models are similar in behavior 
for both u and d quarks, though the agreements in the 
magnitudes are not exact. In AdS/QCD, the nature of 
E(x, bx) for u and d quarks is almost similar when 
plotted against x for fixed values of impact parameter 
b =| bx |, whereas in the diquark model, it shows quite 
different behavior for both u and d quarks. But the 
behaviors of H(x,h±) with respect to x are distinctly 
different for u and d quarks in both the models as can 
be seen in Fig. [8] It is interesting to note that in both 
cases, the GPD i7(x,bx) is larger for u-quarks than 
d-quarks whereas the magnitude of the GPD E(x 1 b±) 
is marginally larger for d-quarks than the that for u- 
quarks at small values of impact parameter b. The sim¬ 
ilar behavior of the GPDs of a phenomenological model 
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Fig. 5 (Color online) Plots of (a) H' u (x,^,t) vs x and —t; (b) the same as in (a) but for d quark; (c) E u {x,(,,t) vs x and —t; 
(d) the same as in (c) but for d quark for fixed £ = 0.15. For £ = 0.15, the minimum value of —t w 0.024 GeV 2 . 


was observed in [35] . Another interesting behavior of all 
the GPDs is that the width of all the distributions in 
transverse impact parameter space decreases as x in¬ 
creases, which implies that the distributions are more 
localized near the center of momentum for higher values 
of x. 

The skewness dependent GPDs in transverse impact 
parameter space for u and d quarks as functions of x 
and b are shown in Fig {TO] for a fixed value of £ = 0.15. 
Similarly, all GPDs as functions of £ and b for a fixed 
value of x = 0.6 are shown in FigJTT] Though there is 
no divergence at x = Q. in the numerical computations, 
the exact value of x = £ has been omitted for techni¬ 
cal reason. At small value of b, H(x,(, bj_) decreases 
for u quark but increases for d quark with increasing £, 
while E(x,(, bj_) decreases with increasing £ for both 
u and d for a fixed value of x. For the fixed values x 
and low £, the peak of u quark is sufficiently large com¬ 
pare to d for H(x , £, bj_) but for E(x , £, bj_), d quark is 
marginally large compare to u quark. Substantial dif¬ 
ferences are observed in both GPDs between u and d 
quarks when the GPDs are plotted against x for fixed 
values of £ and b. H d (x , £, b) seems to be nonzero at 
x = £ in Fig fTiT c). but it is due to the fact that x = £ 


is not included in the plot. It goes to zero as x —> £. It is 
interesting to note that the peaks of all the distributions 
also become broader as £ increases for a fixed value of 
x. This means that the probability of hitting the ac¬ 
tive quark at a larger transverse impact parameter b 
increases as the momentum transfer in the longitudinal 
direction increases. 


3.2 GPDs in longitudinal impact parameter 
space 

The boost invariant longitudinal impact parameter is 
defined as a = | b~P + which is conjugate to the skew¬ 
ness £, the measure of longitudinal momentum transfer. 
The parameter a was first Introduced in El and it was 
shown that the DVCS amplitude in a QED model of 
a dressed electron shows an interesting diffraction pat¬ 
tern in the longitudinal impact parameter space. Since 
Lorentz boosts are kinematical in the light front, the 
correlation defined in the three dimensional position 
space b± and a is frame independent. It was shown 
in the same simple relativistic spin half system of an 
electron dressed with a photon that the GPDs also ex- 
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(c) (d) 

Fig. 6 (Color online) Plots of (a) H u (x,£,t) vs £ ; (b) the s 
as in (c) but for d quark for fixed x = 0.8. 

hibit the similar diffraction pattern in the longitudinal 
impact parameter space m- Similar diffraction pattern 
was also observed in a phenomenological model for pro¬ 
ton GPDsjTS]. So, it is very interesting to investigate if 
the similar pattern is also observed in this light front 
quark model. The GPDs in longitudinal position space 
are defined as: 

H(x,a,t) = j d(e <p+b /2 H(x,C,t), 

= ^ dCe*" H (xX,t), 

E(x,a,t) = J dC,e <p+b /2 E(x,(,t), 

= 7^f o Cf d(e i ^E(x,Ct). ( 28 ) 

Since we are considering the region ( < x < 1, the 
upper limit of C integration Q is given by (max if % is 
larger than <C max , otherwise by x if x is smaller than 
Cmax where the maximum value of £ for a fixed —t is 
given by 

Cmax = 2M|(\/ 1+ Hi -1 )' (29) 




as in (a) but for d quark; (c) E u (x,£,t) vs £; (d) the same 


In FigOZIl we show the GPDs in longitudinal position 
space cr considering the DGLAP region. We observe 
that the GPDs show diffraction pattern in longitudi¬ 
nal impact parameter space, similar to the nature of 
a dressed electron in QED or in a holographic model 
for the meson m- This effect has also been observed 
for the GPDs of a phenomenological model m as well 
as for the chiral odd GPDs of light-front QED model 
[12] . Except for H d (x, cr,t), all the distributions have 
a primary maximum at a = 0 followed by a series of 
secondary maxima. H d (x,a,t) has a peculiar behavior 
having a maximum at cr = 0 for very small —t and 
shows diffraction pattern while for relatively larger val¬ 
ues of —t it shows a minimum at a = 0. The minima in 
E(x, cr, t) occur at the same positions for both u and d 
quarks. In all cases, the position of the first minimum 
moves in to smaller values of cr as —t increases. The 
characteristics of E(x,a,t) for both u and d quarks is 
almost same whereas for H(x, cr, f), the nature of u and 
d quark changes as —t increases. In El. similar diffrac¬ 
tion patterns were observed for DVCS amplitudes with 
both 2 —> 2 and 3—^1 contributions, so we expect 
that the pattern will survive if higher Fock sectors are 
included in the model. 
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Fig. 7 (Color online) Plots of (a) H u (x,(,t) vs x (b) the same as in (a) but for d quark (c) E u (x,C,,t) vs x (d) the same as 
in (c) but for d quark for fixed —t = 0.8 GeV 2 and different values of £• 


4 Transverse charge and magnetization 
densities 


The two dimensional Fourier transform of the Dirac 
form factor gives the transverse charge density in the 
transverse plane for the unpolarized nucleons, 

- / §^ F ^ 2 y ,±M± 

= r ^QMQh) F i(Q 2 ), (30) 

where b represents the impact parameter and Jo is the 
cylindrical Bessel function of order zero. We can write a 
similar formula for charge density for flavor p q ch (b) with 
Fi is replaced by F q . In a similar fashion, one defines 
the magnetization density in the transverse plane by 
the Fourier transform of the Pauli form factor, 

= C ^QMQV)F 2 (Q 2 ), (31) 

whereas, 


Pm(b) = = bj™ ^Q 2 Ji(Qb)F 2 (Q 2 ), (32) 


has the interpretation of anomalous magnetization den¬ 
sity [46] . Since these quantities are not directly mea¬ 
sured in experiments, actual experimental data are not 
available. In Ref. m , an estimation of the proton charge 
and magnetization densities has been done from exper¬ 
imental data of electromagnetic form factors. To get an 
insight into the contributions of the different flavors, 
we evaluate the charge and anomalous magnetization 
densities for the u and d quarks. 

We can define the decompositions of the transverse 
charge and magnetization densities for nucleons in the 
similar way as electromagnetic form factors. The charge 
densities decompositions in terms of two flavors can be 
written as 


P*ch — e uPfch + e dPfchi 

Pch = e uP d fch + e dP 1 fchi (33) 


where e u and are charge of u and d quarks respec¬ 
tively. Due to the charge and isospin symmetry the u 
and d quark densities in the proton are the same as 
the d and u densities in the neutron gang. Under the 
charge and isospin symmetry, we can write 


Pch(b) = p P ch + 


Pch 

2 


Pfch 

2 


db 
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(a) x (b) b [GeV 1 ] 




(c) x (d) b [GeV" 1 ] 

Fig. 8 (Color online) Plots of (a) H u (x, bj_) vs x for fixed values of impact parameter b =| bj_ |; (b) H u (x, bj^) vs b for fixed 
x\ (c) and (d) are the same as in (a) and (b) but for d quarks, b is in GeV -1 . 


P d ch(b)=P P ch + 2p n ch = p d fch , (34) 

where p q ch (b) is the charge density of each quark and 
P}ch i s charge density for each flavor. We can sim¬ 
ilarly decompose p m into magnetization densities for 
each flavor. The flavor contributions to proton charge 
and magnetization densities are e u /dp U / c h an( I e u/dP/m ■ 
Similarly for neutron, the flavor contributions are (e^/u 



In Fig. 1131 we show the charge and anomalous mag¬ 
netization densities for proton and neutron. The plots 
suggest that the light-front diquark model’s results for 
the charge and magnetization density of proton and the 
magnetization density of neutron are in excellent agree¬ 
ment with the two different global parameterizations of 
Kelly QjB] and Bradford at el [T7]. Though the diquark 
model is unable to reproduce the data for the neutron 
charge density at small b, still it is better than the 
AdS/QCD Model-I predictions presented in Ref. [15] . 
In Fig lblT c). one can notice a negatively charged core 
surrounded by a ring of positive charge density for neu¬ 
tron. In Figs. ITH a) and fCT b). we show the charge and 
anomalous magnetization densities for u quark. Simi¬ 
larly for the d quark, the transverse densities are shown 


in Figs. □He) and lldl d). The charge density for d quark 
in diquark model deviates at small b form the two global 
parameterizations of Kelly and Bradford but is in ex¬ 
cellent agreement for u quark. Again, diquark model 
provides better result than the AdS/QCD Model-I re¬ 
sults presented in Ref. m- The anomalous magnetiza¬ 
tion densities in both u and d quarks in the LF diquark 
model match very well with the parameterizations 
For transversely polarized nucleon, the charge den¬ 
sity in the transverse plane is given by m 

Pt{V) = Pch - sin(</>t, - </>s)^T Pm, (35) 

where M is the mass of nucleon and the transverse po¬ 
larization of the nucleon is given by S± = ( cos<f> s x + 
sin (j) s y) and the transverse impact parameter b± = 
b{ cos (f>bX + sin fay). Without loss of generality, the po¬ 
larization of the nucleon is taken along x-axis ie., <p s = 
0. The second term in Eo. (l35l) . provides the deviation 
from circular symmetry of the unpolarized charge den¬ 
sity [47]. The charge densities for the transversely polar¬ 
ized proton and neutron have been shown in Fig. I15f a) 
and M b) . We show the u and d quark charge densities 
for the transversely polarized nucleon in Fig llSI c) and 
fl5ld). Again, the densities in the LF diquak model are 
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(a) 





(b) 


b [GeV~ 


(c) x (d) b [GeV 1 ] 

Fig. 9 (Color online) Plots of (a) £J“(a:,bj_) vs x for fixed values of b =| bj_ |; (b) E u (x, bi) 
the same as in (a) but for d-quarks; and (d) the same as in (b) but for d-quarks. b is in GeV -1 



^ vs h fnr fi-yprl valiifis 




in good agreement with the global parameterizations. 
The comparison of charge densities for the transversely 
polarized and unpolarized proton is shown in Fig. m 
a) and (b) and the similar plots for neutron are shown 
in FigHni(c) and (d). For the nucleons polarized along 
the x direction, the charge densities get shifted towards 
negative b y direction for proton. The deviation is much 
larger for the neutron compared to the proton. Due to 
large anomalous magnetic moment which produces an 
induced electric dipole moment in y-direction, the dis¬ 
tortion shows a dipolar pattern in the case of neutron 
m- The behaviors are in agreement with the results 
reported in Refs. mmm- 


We compare the up quark charge densities for the 
transversely polarized and unpolarized nucleon in Fig llTI 
(a) and (b) and the similar plots for d quark are shown 
in Fig. E] (c) and (d). The deviation or distortion from 
the symmetric unpolarized density is more for down 
quark than the up quark. For the nucleons polarized in 
x direction, the charge density shifts towards positive 
by direction for d quark but in opposite direction for 
the u quark. 


5 Summary and conclusions 

The parameters in a light front quark-diquark model 
of nucleons [I] are found to be inconsistent with the 
experimental data. We have re-evaluated the parame¬ 
ters in this model for the AdS/QCD scale parameter 
k = 0.4066 GeV which was previously obtained by fit¬ 
ting the nucleon form factors in soft wall AdS/QCD [HJ 
The new parameters reproduce the experimental 
data for the nucleon form factor quite well for a wide 
range of Q 2 values. We have compared our results with 
AdS/QCD soft wall model. Then, we have evaluated 
the GPDs for u and d quark in proton for both zero 
and nonzero skewness in the light front quark-diquark 
model. We observed that all the GPDs in the momen¬ 
tum space as well as in the transverse impact parameter 
space are more or less in agreement with the results of 
AdS/QCD. We have calculated the GPDs for nonzero 
skewness in the DGLAP region(i.e., for x > </). The 
peaks of the distributions move to higher values of x 
for fixed C with increasing —t. In the model, the behav¬ 
iors of the GPD H in impact parameter space for u and 
d quarks are quite different when plotted in both x and 
b. The difference in the behaviors of E(x, t ) for u and 
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Fig. 10 (Color online) Plots of (a) H u (x,(,h±) vs x and b =| h± |; (b) E u (x,(,h±) vs. x and 6; (c) the same as in (a) but 
for d-quarks; and (d) the same as in (b) but for d-quarks for fixed value of £ = 0.15. 


d quarks are clearly observed when plotted against x. 
For nonzero skewness, we have also shown the GPDs 
in longitudinal impact parameter space a. We found 
that both the GPDs H and E for u and d quarks in 
a space show diffraction patterns. Similar diffraction 
patterns also have been observed in some other mod¬ 
els. In case of E(x,a,t), the qualitative nature of the 
diffraction pattern is same for both u and d quarks. For 
H d {x, a,t), the diffraction pattern is observed only for 
small —t values; as —t increases, a dip appears at the 
centre( at a = 0). 

Finally, we have presented the transverse charge and 
magnetization densities for nucleon and also for indi¬ 
vidual quarks. The results are consistent with two phe¬ 
nomenological narameterizations [161fT7] . The unpolar¬ 
ized densities are axially symmetric whereas the charge 
densities get distorted for transversely polarized nu¬ 
cleon. The charge density is shifted along y direction 
if the nucleon is polarized along x direction. The charge 
density for transversely polarized neutron shows a dipole 
pattern. The shift of charge density of u quark for trans¬ 
versely polarized nucleon from the symmetric unpolar¬ 
ized density is smaller than d quark and in opposite 
direction. 


Appendix A: comparison of GPDs in 
quark-diquark model with a double 
distribution(DD) model 


The GPDs for £ = 0 admit a density interpretation 
when one takes the Fourier transform to the impact pa¬ 
rameter space but in experiments, £ is always nonzero. 

In recent past, there have been a lot of works to model 
GPDs with nonzero skewness by modeling relevant DDs 
mm- In this section, we compare our results for nonzero 
skewness with the GPDs modeled from the Double Dis¬ 
tributions (dd) [H21ISnilS3] - The GPDs have an integral 
representation in terms of the double distributions /(/?, a, t). 
For the valance quarks, the GPDs can be written as 
1 j. 1-/3 

d/3 da 5(x-/3-(a) /®(/3,a,f), 

J13-1 



(A.l) 


where F® = HE%. Here, we use the factorized DD 
ansatz for the GPDs as suggested by Musatov and Radyushkin 

E2 

/«(/3,M) = F*(P,0,t)h(p,a), (A.2) 

where the weight function h(f3,a) generates the skew¬ 
ness dependence of the GPDs and satisfies the normal- 
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Fig. 11 (Color online) Plots of (a) H u (x, £, bj_) vs £ and b =| bj_ |; (b) £“(i, £, bx) vs. £ and 6; (c) the same as in (a) but 
for d-quarks; and (d) the same as in (b) but for d-quarks for fixed value of x = 0.6. 


ization condition 
/•1-I/3I 

/ h(/3,a)da = 1. (A. 3) 

J- i+l/3| 

The general form of the profile function is given by |55j 

r(2N + 2) [(1 - \P\) 2 - a 2 ] N 


h (N) (P,a) = 


,(A.4) 


2 2J v+ir 2 (iV+ 1 ) (1 ^ |/3|) 2Ar+1 

where the parameter N governs the width of the func¬ 
tion. We use the N = 2 profile function. The simi¬ 
lar profile function for N = 2 has been used in many 
phenomenological model of DVCS and exclusive meson 
production mrnmm- Inserting Eq. IA.2I in the Eq. 
IAJ1 with the help of delta-function one can perform 
the integral over a and obtains 

3 d(3 

4^3 


= 77s /„ i_p 


W0,i) 


< A - 5 » 


-P 

for x > C, the integration boundaries are 

c 


Pmin — % 


= X + 


1-C 

c 

1-C 


(1 — x ) 
(1 — x). 


(A.6) 


In Fig. [lH] and Fig. [T9l we show the skewness de¬ 
pendent GPDs calculated using double distribution pa¬ 
rameterization and compare with the results directly 
calculated in the quark-diquark model. FigfTBl suggests 
that for small £ and large —t, the results of double dis¬ 
tribution are more or less in agreement with the diquark 
model results, while Figflbl shows that at moderate or 
high values of skewness (, only at higher x the two 
models agree but otherwise the agreement is lost. 
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Fig. 13 (Color online) Plots of transverse charge and anomalous magnetization densities for nucleon, (a) and (b) represent 
p c h and p m for the proton, (c) and (d) same as proton but for neutron. Dashed lines represent the parameterization of Kelly 
[16| . and the lines with circles represent the parameterization of Bradford at el m, dot-dashed lines are for soft-wall model 
m ■ The solid lines represent the light-front scalar diquark model. 
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Fig. 14 (Color online) Plots of transverse charge and anomalous magnetization densities for quark, (a) and (b) represent p c h 
and p m for the u quark, (c) and (d) same as u but for d quark. Dashed lines represent the parameterization of Kelly [16] . and 
the lines with circles represent the parameterization of Bradford at el m, dot-dashed lines are soft-wall model Il8l . The solid 
lines represent the light-front scalar diquark model. 
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Fig. 15 (Color online) The charge densities for the transversely polarized (a) proton (b) neutron and (c) up (d) down quark 
charge densities for the transversely polarized nucleon. Dashed line represents the parameterization of Kelly m, and line 
with circles represents the parameterization of Bradford at el HU, dot-dashed line is soft-wall Model-I in [18] . The solid line 
represents this work for light-front scalar diquark model. 
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(c) (d) 

Fig. 16 (Color online) The charge densities in the transverse plane for the (a) unpolarized proton (b) transversely polarized 
proton and (c) unpolarized neutron and (d) transversely polarized neutron. Transverse polarization is along x direction. 
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(c) (d) 

Fig. 17 (Color online) The charge densities in the transverse plane of u quark for (a) unpolarized (b) transversely polarized 
nucleon and d quark for (c) unpolarized (d) transversely polarized nucleon. Polarization is along ^-direction. 
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(a) X (b) 




Fig. 18 (Color online) Plots of (a) H u (x, C, t) vs x (b) the same as in (a) but for d quark (c) E u (x, C, t) vs x (d) the same as 
in (c) but for d quark for fixed ( — 0.15 and three different values of —t = 0.3GeV 2 , l.OGeV 2 and 3.0GeV 2 . Red dashed lines 
represent the results calculated using double distribution for the same values of £ and — t. 





















22 


(a) 


(c) 



Fig. 19 (Color online) Plots of (a) H u {x,C,,t) vs x (b) the same as in (a) but for d quark (c) E u {x,C,,t ) vs x (d) the same as 
in (c) but for d quark for fixed — t = 1.0 GeV 2 and three different values of C = 0.05, 0.2 , 0.4. Red dashed lines represent the 
results calculated using double distribution for same values of — t and £(. 









